Volatility Effects on the Escape Time in Financial Market Models 
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We shortly review the statistical properties of the escape times, or hitting times, for stock 
price returns by using different models which describe the stock market evolution. We 
compare the probability function (PF) of these escape times with that obtained from real 
market data. Afterwards we analyze in detail the effect both of noise and different initial 
conditions on the escape time in a market model with stochastic volatility and a cubic 
nonlinearity. For this model we compare the PF of the stock price returns, the PF of 
the volatility and the return correlation with the same statistical characteristics obtained 
from real market data. 
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I. INTRODUCTION 

Econophysics is a developing interdisciplinary 
research field in recent years. It applies theo- 
ries and methods originally developed by physi- 
cists in statistical physics and complexity in or- 
der to solve problems in economics such as those 
strictly related to the analysis of financial market 
data [Anderson et al, 1988, 1997; Mantegna & 
Stanley, 2000; Bouchaud & Potters, 2004]. Most 
of the work in econophysics has been focused on 
empirical studies of different phenomena to dis- 
cover some universal laws. Recently more effort 
has been done to construct new models. In a real 
market the stock option evolution is determined 
by many traders which interact with each other 
and use different strategy to increase their own 
profit. The market is then "pushed" by many 
different forces, which often affect the system 
in such a way that every deterministic forecast 
is impossible. In fact people act in the market 
so that any forecast results to be unpredictable. 
The arbitrariness of each choice, together with 
the non-linearity of the system, leads to consider 
the stock option market as a complex system 
where the randomness of the human behavior 
can be modelled by using stochastic processes. 

For decades the geometric Brownian motion. 
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proposed by Black and Scholes [Black & Scholes, 
1973] to address quantitatively the problem of 
option prices, was widely accepted as one of 
the most universal models for speculative mar- 
kets. However, it is not adequate to correctly 
describe financial market behavior [Mantegna &; 
Stanley, 2000; Bouchaud & Potters, 2004]. A 
correction to Black-Scholes model has been pro- 
posed by introducing stochastic volatility mod- 
els. These models are used in the field of quan- 
titative finance to evaluate derivative securities, 
such as options, and are based on a category of 
stochastic processes that have stochastic second 
moments. In finance, two categories of stochas- 
tic processes are widely used to model stochas- 
tic second moments. One is represented by the 
stochastic volatility models (SVMs), the other 
one by ARCH/GARCH models [Engle, 1982; 
Bollerslev, 1986], where the present volatility 
depends on the past values of the square re- 
turn (ARCH) and also on the past values of the 
volatility (GARCH). Both ARCH/GARCH and 
stochastic volatility models derive their random- 
ness from white noise processes. The difference 
is that an ARCH/GARCH process depends on 
just one white noise, while SVMs generally de- 
pend on two white noises and they model the 
tendency of volatility to revert to some long- 
run mean value. Stochastic volatility models 
address many of the short-comings of popular 



option pricing models such as the Black-Scholes 
model [Black & Scholes, 1973] and the Cox- 
Ingersoll-Ross (CIR) model [Cox et al, 1985]. 
In particular, these models assume that the un- 
derlier volatility is constant over the life of the 
derivative, and unaffected by the changes in the 
price level of the under lier. However, these 
models can not explain long-observed anoma- 
lies such as the volatility smile [Fouque et al., 
2000] and some stylized facts observed in finan- 
cial time series such as long range memory and 
clustering of the volatility, which indicate that 
volatility does tend to vary over the time [Da- 
corogna et al, 2001]. By assuming that the 
volatility of the underlying price is a stochas- 
tic process rather than a constant parameter, 
it becomes possible to more accurately model 
derivatives. In SVMs the volatility is changing 
randomly according to some stochastic differen- 
tial equation or some discrete random processes. 
Recently, models of financial markets reproduc- 
ing the most prominent statistical properties of 
stock market data, whose dynamics is governed 
by non-linear stochastic differential equations, 
have been proposed [Malcai et al., 2002; Bor- 
land, 2002; Borland, 2002b; Hatchett & Kiihn, 
2006; Bouchaud & Cont, 1998; Bouchaud, 2001; 
Bouchaud, 2002; Sornette, 2003; Bonanno et al., 
2006; Bonanno et al., 2007]. 

In particular some models have been used 
where the market dynamics is governed, close 
to crisis period, by a cubic potential with a 
metastable state [Bouchaud & Potters, 2004; 
Bouchaud & Cont, 1998; Bouchaud, 2001; 
Bouchaud, 2002; Bonanno et al., 2006; Bonanno 
et al., 2007]. The metastable state is connected 
with the stability of normal days, when the fi- 
nancial market shows a normal behavior. Con- 
versely the presence of a crisis is modelled as 
an escape event from the metastable state and 
the subsequent trajectory. The importance of 
metastable states in real systems, ranging from 
biology, chemistry, ecology to population dy- 
namics, social sciences, economics, caused re- 
searchers to devote many efforts to investigate 
the dynamics of metastable systems, finding that 
they can be stabilized by the presence of suitable 
levels of noise intensity [Mantegna & Spagnolo, 
1996; Mielke, 2000; Agudov & Spagnolo, 2001; 



Dubkov et al, 2004; Fiasconaro et al, 2005; Fi- 
asconaro et al, 2006; Bonanno et al, 2007]. 

Our focus in this paper is to analyze the sta- 
tistical properties of the escape times in differ- 
ent market models, by comparing the probability 
function (PF) with that observed in real market 
data. Recent work has been done on the mean 
exit time [Bonanno & Spagnolo, 2005; Montero 
et al, 2005] and the waiting time distribution 
in financial time series [Raberto et al., 2002]. 
Here, starting from the geometric random walk 
model we shortly review the statistical proper- 
ties of the escape times for stock price returns 
in some stochastic volatility models as GARCH, 
Heston and nonlinear Heston models. In the last 
model, recently proposed by the authors [Bo- 
nanno et al, 2006; Bonanno et al, 2007] and 
characterized by a cubic nonlinearity, we com- 
pare some of the main statistical characteristics, 
that is the PF of the stock price returns, the 
PF of the volatility and the return correlation, 
with the same quantities obtained from real mar- 
ket data. We also analyze in detail the effect 
of the noise and different initial conditions on 
the escape times in this nonlinear Heston model 
(NLH). 



II. ESCAPE TIMES IN STOCK MARKET 
MODELS 

The average escape time of a Brownian parti- 
cle, moving in a potential profile, is a well-known 
problem in Physics [Gardiner, 2004]. This quan- 
tity is defined as < T > = f{t)dt, where f{t) 
is the probability function of the escape events 
from a certain region of the potential (see Fig.[l]). 
In other words, the average escape time is the 
mean time that a particle, starting from a cer- 
tain initial point, takes to pass through a given 
threshold. This average time, obtained by us- 
ing the backward Fokker-Planck equation, corre- 
sponds to the first passage time used in statisti- 
cal physics [Redner, 2001; Inoue & Sazuka, 2007] 
and to the first hitting time defined in econo- 
physics [Bonanno & Spagnolo, 2005; Montero et 
al., 2005]. 
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FIG. 1 Cubic potential used in the dynamical equa- 
tion for the process x(t) of Eq. (7). The white circles 
in the figure indicate the starting positions used in 
our simulations. 

A. The geometric Brownian motion model 

A common starting point for many theories in 
economics and finance is that the stock price, in 
the continuous limit, is a stochastic multiplica- 
tive process defined, in the Ito sense, as 

dp{t)='yp{t)-dt + a-p{t)-dW{t) (1) 

where 7 and a are, in the market dynamics, the 
expected average growth for the price and the 
expected noise intensity {volatility) respectively. 
The price return dp/p = dlnp{t) obeys the fol- 
lowing additive stochastic differential equation 
2 

dlnp{t) = {ii-^)-dt + a-dW{t) . (2) 

This simple market model, proposed by Black 
and Scholes, catches one of the most important 
stylized facts of the financial markets, that is the 
short range correlation for price returns. This 
characteristic is necessary in order to warrant 
market efficiency. 

In the geometric Brownian motion the re- 
turns are independent on each other, so the prob- 
ability to observe a value after a certain barrier 
is given by the probability that the "particle" 
doesn't escape after n — 1 time steps, multiplied 
by the escape probability at the n*^ step 

F{t) = (3) 
= (1 — p) • exp [{n — 1) Inp] , n = — 



where p is the probability to observe a return in- 
side the region limited by the barrier, Ai is the 
observation time step and r is the escape time. 
So the behavior of the PF of hitting times is ex- 
ponential. The geometric Brownian motion how- 
ever is not adequate to describe financial market 
behavior. In fact the volatility is considered as 
a constant parameter and the PF of the price 
is a log-normal distribution. As a consequence 
many observations of real data are in clear dis- 
agreement with this model [Mantegna & Stanley, 
2000; Bouchaud & Potters, 2004]. 

B. Stocliastic volatility models 

1. The GARCH model 

Data on financial return volatility are influ- 
enced by time dependent information flow which 
results in pronounced temporal volatility cluster- 
ing. These time series can be parameterized us- 
ing Generalized Autoregressive Conditional Het- 
eroskedastic (GARCH) models. It has been 
found that GARCH models can provide good 
in-sample parameter estimates and, when the 
appropriate volatility measure is used, reliable 
out-of-sample volatility forecasts [Anderson & 
Bollerslev, 1998]. 

The GARCH(p,q) process, which is essen- 
tially a random multiplicative process, is the 
generalization of the ARCH process and com- 
bines linearly the present return with the p previ- 
ous values of the variance and the q previous val- 
ues of the square return [Bollerslev, 1986]. The 
process is described by the equations 

a^ = ao + aiXt_i H h OgX^-g + /3i<7f-i 

+ \- PpO-t-p, xt = r]fat, (4) 

where Oj and Pi are parameters that can be es- 
timated by means of a best fit of real market 
data, r]t is an independent identically distributed 
random process with zero mean and unit vari- 
ance. Using the assumption of Guassian con- 
ditional PF, rjt is Gaussian. In Eq.([4]) xt is a 
stochastic process representing price returns and 
is characterized by a standard deviation at- The 
GARCH process has a non-constant conditional 
variance but the variance observed in long time 
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period, called unconditional variance, is instead 
constant and can be calculated function of 
the model parameters. We shall consider the 
simpler GARCH(1,1) model 

CTt = ao + (aiVt-i + /3i)<7?-i^ xt = r]f at . (5) 

The autocorrelation function of the process xt 
is proportional to a delta function, while the 
process x^ has a correlation characteristic time 
equal to r =[ ln{ai + l""*^ and the uncondi- 
tional variance equal to = ao/(l — ai — 
By a fitting procedure between the previous ex- 
pressions of r and and the empirical values of 
the same quantities, we can easily estimate the 
three parameters ao, ai and /?i which character- 
ize the model. 

2. The Heston model 

The Heston model introduced by Heston [He- 
ston, 1993] is a commonly used stochastic volatil- 
ity model. It received a great attention in the 
financial literature specially in connection with 
option pricing [Fouque et al, 2000]. The He- 
ston model was verified empirically with both 
stocks [Silva & Yakovenko, 2003, Dragulescu & 
Yakovenko, 2002] and options [Huh & White, 
1987; Hull, 2004], and good agreement with the 
data has been found. It was also recently in- 
vestigated by econophysicists [Micciche et a/., 
2002; Dragulescu & Yakovenko, 2002; Silva et 
ai, 2004; Bonanno & Spagnolo, 2005]. The 
model is defined by two coupled stochastic dif- 
ferential equations which represent the stock dy- 
namics by the log-normal geometric Brownian 
motion stock process and the Cox-Ingersoll-Ross 
(CIR) mean-reverting process (SDE), first intro- 
duced to model the short term interest rate [Cox 
et a/., 1985]. By considering the log of the price 
x{t) = In p{t) the SDEs of the model are 



dx{t) = i^n - ^ j ■ dt + ■ dWi{t) 

dv{t) = a{b - v{t)) ■ dt + c^/v{t) ■ dWc{t) 
dW^{t) = p-dWi{t) + ^/l- ■dW2{t), (6) 

where p is the trend of the market, Wi{t) and 
W2{t) are uncorrelated Wiener processes with 



the usual statistical properties < dWi > = 0, 
< dWi{t)dWj{t') >= dt6{t-t') 6ij = 1,2), 
and p is the cross correlation coefficient between 
the noise sources. Here v is the CIR process, 
which is defined by three parameters: b, a and c. 
They represent respectively the long term vari- 
ance, the rate of mean reversion to the long term 
variance, and the volatility of variance, often 
called the volatility of volatility. The stochas- 
tic volatility v{t) is characterized by exponential 
autocorrelation and volatility clustering [Cont, 
2001; Bouchaud & Potters, 2004; Bonanno et al, 
2006], that is alternating calm with burst peri- 
ods of volatility. 

We end this paragraph comparing the PF of 
the escape times r of the price returns obtained 
from real market data with the PFs obtained by 
using the three previous models. We use a set 
of returns obtained from the daily closure prices 
for 1071 stocks traded at the NYSE and continu- 
ously present in the 12— year period 1987— 1998 
(3030 trading days). From this data set we ob- 
tain the time series of the returns and we cal- 
culate the time to hit a fixed threshold starting 
from a fixed initial position. The parameters in 
the models were chosen by means of a best fit, 
in order to reproduce the correlation properties 
and the variance appropriate for real market [Bo- 
nanno et al, 2005]. We choose two thresholds 
to define the start and the end for the random 
walk. Specifically we calculate the standard de- 
viation cjj, with i = 1, . . . , 1071 for each stock 
over the whole 12-year period. Then we set the 
initial threshold at the value —0.1 • cjj and as fi- 
nal threshold the value —2 • ai. The thresholds 
are different for each stock, the final threshold is 
considered as an absorbing barrier. The results 
of our simulations together with the real market 
data are shown in the following Fig. [2l In this fig- 
ure the exponential behavior represents the PF 
of the escape times for the geometric Brownian 
motion model, which is not adequate do describe 
correctly the PF of r over the entire time axis. 
The GARCH model provides a better qualitative 
agreement with real data for lower escape times 
and gives the exponential behaviour in the re- 
gion of large escape times. Here the geometric 
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Brownian model reproduces well the real data^, 
whereas the Heston model is able to reproduce 
almost entirely the empirical PF. 
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FIG. 2 Probability function (PF) of escape times of 
stock price returns for: real market data (circle), ge- 
ometric Brownian motion model (black solid line), 
GARCH model (black broken line) and Heston model 
(diamond) . 



C. The nonlinear Heston model (NLH) 

To consider feedback effects on the price fluc- 
tuations and different dynamical regimes, simi- 
larly to what happens in financial markets dur- 
ing normal activity and in special days with rel- 
atively strong variations of the price [Bouchaud 
& Potters, 2004; Bouchaud & Cont, 1998; 
Bouchaud, 2001; Bouchaud, 2002; Bonanno et 
al, 2006; Bonanno et al, 2007], we proposed 
a generalization of the Heston model, by con- 
sidering a cubic nonlinearity in the SDE of the 
log of the price x{t) = In p{t) (first equation 
in (6)) [Bonanno et al, 2006; Bonanno et al, 
2007]. This nonlinearity allows us to describe 
these different dynamical regimes by the motion 
of a fictitious " Brownian particle''^ moving in an 
effective potential with a metastable state. The 



equations of the new model are obtained by re- 
placing in Eqs. ([6]) the parameter // with the neg- 
ative derivative of the nonlinear cubic potential 



dx{t) 

dv{t) 
dW,{t) 



^ + ^) dt + ^)dW,{t){l) 



a{h - v{t))dt + cy/v{t)dW^{t) 
pdWi{t) + ^l- p'^dW2{t), 



(8) 
(9) 



^ By changing the fit parameters a\ and /3i for the 
GARCH model it is possible to obtain a better good 
agreement with the real data. 



where U{x) = 2x^ + 3x^ is the effective cubic 
potential with a metastable state at Xme = 0, 
a maximum at Xm = — 1, and a cross point be- 
tween the potential and the 
(see Fig. [1]). The average exit time of the sys- 
tem from the stable to the unstable domain of 
the potential shown in Fig. [1] may be prolonged 
by imposing external noise: this phenomenon is 
named noise enhanced stability (NES). The sta- 
bility of systems with a metastable state can 
be increased by enhancing the lifetime of the 
metastable state or the average exit time of the 
system from the well. The NES effect was ex- 
perimentally observed in a tunnel diode [Man- 
tegna & Spagnolo, 1996] and in an underdamped 
Josephson junction [Sun et el., 2007] and theo- 
retically predicted in a wide variety of systems 
such as chaotic map, Josephson junctions, neu- 
ronal dynamics models and tumor-immune sys- 
tem models [Mielke, 2000; Agudov & Spagnolo, 
2001; Dubkov et al, 2004; Pankratova et al, 
2004; Pankratov & Spagnolo, 2004; Fiasconaro 
et al, 2005; Fiasconaro et al, 2006]. Two dif- 
ferent dynamical regimes are observed depend- 
ing on the initial position of the Brownian par- 
ticle along the potential profile. One is charac- 
terized by a nonmonotonic behavior of the life- 
time, as a function of the noise intensity (here 
the volatility v{t)), for initial positions Xq < Xc- 
The other one features a divergence of the life- 
time when the noise tends to zero for initial posi- 
tions Xc < Xo < Xm , implying that the Brownian 
particle remains trapped inside the metastable 
state in the limit of small noise intensities. In 
this dynamical regime a nonmonotonic behav- 
ior of the lifetime with a minimum and a maxi- 
mum as a function of the noise intensity is also 
observed. This trapping phenomenon is always 
observable when initial unstable positions of the 
Brownian particle are near a metatastable state 
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of the system investigated. The NES effect and 
its different dynamical regimes can be explained 
considering the barrier "seen" by the Brownian 
particle starting at the initial position Xq, that 
is AUin = U{xmax) - U{xo), and by compar- 
ing it with the height of the barrier AU charac- 
terizing the metastable state (see Fig. [TJ [Agu- 
dov & Spagnolo, 2001; Fiasconaro et al, 2005]. 
For example for unstable initial positions such 
as Xc < Xo < Xm we have AUm < AU and from 
a probabilistic point of view, it is easier to en- 
ter into the well than to escape from, once the 
particle is entered. So a small amount of noise 
can increase the lifetime of the metastable state. 
When the noise intensity v is much greater than 
AU, the typical exponential behavior is recov- 
ered. 

By investigating the mean escape time 
(MET), as a function of the model parameters a, 
b and c, we found the parameter region where a 
nonmonotonic behavior of MET is observable in 
our NLH model with stochastic volatility v{t) 
[Bonanno et al, 2006; Bonanno et al, 2007]. 
This behaviour is similar to that observed for 
MET versus v in the NES effect with constant 
volatility v. We call the enhancement of the 
mean escape time (MET), with a nonmonotonic 
behavior as a function of the model parameters, 
NES effect in a broad sense. Two limit regimes 
characterize our NLH model, one corresponding 
to the case a = 0, with only the noise term in 
the equation for the volatility v{t), and the other 
one corresponding to the case c = with only the 
reverting term in the same Eq. (8). In the first 
case (a = 0), the system becomes too noisy and 
the NES effect is not observable in the behavior 
of MET clS Si function of the parameter c. In the 
second case (c = 0), after an exponential tran- 
sient, the volatility reaches the asymptotic value 
b, and the NES effect is observable as a function 
of b. This case corresponds to the usual constant 
volatility regime. 

By considering the two noise sources Wi{t) 
and Wc{t) of Eqs. (7) and (8) completely uncor- 
related (/? = 0) , the results of simulations of the 
NLH model (Eqs. (7)-(9)), in the second case 
(c = 0), are reported in Fig. [3l where MET ver- 
sus b is plotted for three different starting un- 
stable initial positions and for c = 0. The sim- 
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FIG. 3 Mean escape time (MET) for 3 different un- 
stable starting positions, when only the reverting 
term is present: a = 10^^, c = 0. The curves are 
averaged over 10^ escape events. 



ulations were performed considering the initial 
positions of the process x{t) in the unstable re- 
gion [xc, Xrri\ and using an absorbing barrier at 
X = —6.0. When the process x{t) hits the bar- 
rier, the escape time is registered and another 
simulation is started, placing the process at the 
same starting position Xq, but using the volatil- 
ity value of the barrier hitting time. The non- 
monotonic behavior, which is more evident for 
starting positions near the maximum of the po- 
tential, is always present. After the maximum, 
when the values of b are much greater than the 
potential barrier height, the exponential behav- 
ior is recovered. The results of our simulations 
show that the NES effect can be observed as a 
function of the volatility reverting level 6, the ef- 
fect being modulated by the parameter (ab)/c. 
The phenomenon disappears if the noise term is 
predominant in comparison with the reverting 
term. Moreover the effect is no more observable 
if the parameter c pushes the system towards a 
too noisy region. When the noise term is cou- 
pled to the reverting term, we observe the NES 
effect as a function of the parameter c. The ef- 
fect disappears if b is so high as to saturate the 
system. 

In financial markets the log of the price x{t) = 
In p{t) and the volatility v{t) can be correlated 
(/9 7^ 0), and a negative correlation between the 
processes is known as leverage effect [Fouque et 
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al., 2000]. A negative correlation between the 
logarithm of the price and the volatility means 
that a decrease in x{t) induces an increase in 
the volatility v{t), and this causes the Brown- 
ian particle to escape easily from the well. As a 
consequence the mean lifetime of the metastable 
state decreases, even if the nonmonotonic behav- 
ior is still observable. On the contrary, when 
the correlation p is positive, decrease in x{t) in- 
deed is associated with decrease in the volatil- 
ity, the Brownian particle therefore stays more 
inside the well. The escape process becomes 
slow and this increases further the lifetime of 
the metastable state. The presence of correla- 
tion between the stochastic volatility and the 
noise source which affects directly the dynamics 
of the quantity x{t) = Inp(t) (as in usual mar- 
ket models) can influence therefore the stability 
of the market. Specifically a positive correla- 
tion between x{t) and volatility v{t) slows down 
the walker escape process, that is it delays the 
crash phenomenon increasing the stability of the 
market. Conversely a negative correlation accel- 
erates the escape process, lowering the stability 
of the system [Bonanno et al., 2006]. 



III. ROLE OF THE INITIAL CONDITIONS AND 
STATISTICAL FEATURES 

In this last section we study, for the uncor- 
related {p = 0) NHL model (Eqs. (7)- (9)), the 
role of the initial position of the fictitious Brow- 
nian particle on the mean escape time (MET). 
In particular we fix an escape barrier (thresh- 
old) and we analyze the behaviour of MET as 
a function of a, b and c (CIR process param- 
eters) for different start positions (see Fig [1]). 
First we consider the behavior of MET func- 
tion of the reverting term b. In Fig. [J] (panel 
(a) of part (A)) we show the curves averaged 
on 10'' escape events. The curves inside all the 
other panels have been obtained averaging on 
10^ realizations. In our simulations we consider 
two different values of the parameter c, namely 
c = 10~^,10, and eight values of the parame- 
ter a, that is: (A) a = 10"^ lO'^, lO""^, 10"^ 
and (B) a = 10-\ 1, 10, 10^. Each panel cor- 
responds to a different value of initial position. 



namely: (a) Xq = —0.75, (b) Xq = —1.10, (c) 
Xo = —1.40, (d) Xo = —1.60. Inside each panel 
different curves correspond to different values of 
a. First of all we comment the panels (b), (c) and 
(d) of part (A), related to unstable initial posi- 
tions outside the potential well. The nonmono- 
tonic shape, characteristic of the NFS effect, is 
clearly shown in these three panels. This behav- 
ior is shifted towards higher values of b as the 
parameter a decreases, and it is always present. 
The NFS effect is more pronounced for initial po- 
sitions near the top of the potential barrier. For 
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FIG. 4 Mean escape time (MET) as a function of 
reverting level b. (A) c = 10~^, (B) c = 10. Each 
panel corresponds to different values of initial posi- 
tion. Inside each panel different curves correspond to 
different values of a. 

initial positions far from the maximum of the 
potential, the trapping event becomes less prob- 
able. To obtain a more pronounced NFS effect 
we should consider very low values of the absorb- 
ing barrier, that is x <C —6, which are meaning- 
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less from financial market point of view. For the 
higher value c = 10, we observe the nonmono- 
tonic behavior only for a very great value of the 
parameter a, that is for a = 100 3> c (see panels 
(b), (c), (d) of Fig. [3J3). For further increase of 
the parameter c the noise experienced by the sys- 
tem is much greater than the effective potential 
barrier " seen" by the fictitious Brownian particle 
and the NFS effect is never observable. We note 
that the parameters a and c play a regulatory 
role in Fq. (8). In fact for a ^ c the drift term 
is predominant while for a <C c the dynamics is 
driven by the noise term, unless the parameter b 
takes great values. The nonmonotonic behavior 
is observed for a <C c, provided that b ^ c. For 
increasing values of a the system approaches the 
revert-only regime and we recover the behavior 
shown in Fig. [3j 

Now we consider the panel of Fig |4]A related 
to the unstable initial position outside the po- 
tential well. For very low values of the parame- 
ter a the nonmonotonic behavior is absent and 
the mean escape time (MET) decreases mono- 
tonically with strong fluctuations. We recover 
the similar behavior obtained in the limit case 
of a = and discussed in [Bonanno et al., 2007]. 
For a <C c, in fact, we can neglect the revert- 
ing term in Eq. (8) and the volatility is propor- 
tional to the square of the Wiener process. The 
dynamics is dominated by the noise term with 
large fluctuations for the MET. This behavior is 
mainly due to the presence of the Ito term in 
Eq. (7) for log returns x{t). The Ito term mod- 
ifles randomly the potential shape of Fig. [1] in 
such a way that the potential barrier disappears 
for greater values of the volatility v{t), produc- 
ing a random enhancement of the escape pro- 
cess. Increasing the value of a these fluctuations 
disappear because the reverting term becomes 
more important. This is the behavior shown for 
a = 10~^. A further increase of a causes the 
revert term to dominate the dynamics with re- 
spect to the noise term. Moreover, because the 
initial unstable position Xq = —0.75 is near the 
maximum of the potential well, we recover the 
divergent dynamical regime characterized by a 
nonmonotonic behavior with a minimum and a 
maximum of MET clS Si function of the noise in- 
tensity, here represented by the parameter b [Fi- 



asconaro et al, 2005]. For very low values of b, 
the flctitious Brownian particle is trapped inside 
the potential well with a divergence of MET in 
the limit 6 — > 0. For increasing b the particle can 
escape more easily, and the MET decreases, as 
long as the noise intensity, represented by the pa- 
rameter b, reaches the value 0.15 corresponding 
to the barrier height Af7j„ "seen" by the parti- 
cle. Close to this value of b, the escape process 
is slowed down, because the probability of reen- 
tering the well is equal to that of escaping from. 
This behavior is represented by the minimum of 
MET at 6 ~ 0.15 in the panel (a) of Fig. HA.. 
By increasing b, the particle escaped from the 
well can reenter as long as the noise intensity is 
comparable with the height of the potential bar- 
rier. The MET therefore increases until reach 
a maximum at 6 ~ AJ7 = 1. At higher val- 
ues of b, one recovers a monotonic decreasing 
behavior of the MET. The same nonmonotonic 
behavior with a minimum and a maximum is vis- 
ible in Fig. [3] for Xq = —1.10. Now we consider 

H h a = lO"^ a = 2.7 x lO"^ O — 0= a = 7.3 s lO"^ A A a = <>.<; 




ttr" lo*^ lu** lo' la^ lu" nr'^ iir' lo* m' lo'' icj" 




lu"'^ lo"^ III** lo' HI* lu" la ' 11)"' lo* lu' lo* i[i° 



FIG. 5 Mean escape time (MET) as a function of 
the noise intensity c for a fixed value of 6 (6 = 10~^). 
Each panel corresponds to different values of initial 
condition as in Fig H) Inside each panel different 
curves correspond to the following values of a: cross 
10-3, star 2.7 x lO'^, diamond 7.3 x lO'S triangle 
6.6. 

the dependence of MET on the noise intensity 
c. Fig. [5] shows the curves of MET versus c, av- 
eraged over 10^ escape events in panels (b), (c), 
(d) and 10^ escape events in panel (a). First 
of all we consider the unstable initial positions 
outside the well, that is panels (b), (c) and (d). 
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Each panel corresponds to a different value of 
initial position as in Fig [H Inside each panel 
different curves correspond to different values of 
a. The shape of the curves is similar to that ob- 
served in Fig. m Specifically for small values of 
a, when the reverting term is negligible, the ab- 
sence of the nonmonotonic behavior is expected. 
By increasing a the nonmonotonic behavior is 
recovered. Again the NES effect is more pro- 
nounced for initial positions near the maximum 
of the potential. For initial position inside the 
potential well as in panel (a) , we observe a diver- 
gent behavior of MET for three values of the pa- 
rameter a (namely a = 2.7 • 10~^, 7.3 • 10~^, 6.6), 
because of the small value of 6 = 10~^. Recall 
that the volatility v{t) reverts towards a long 
term mean squared volatility b with relaxation 
time given by a~^. So, for increasing values of a 
the Brownian fictitious particle experiences the 
low value of the noise intensity b = 10~^ in a 
shorter time and therefore the particle is trapped 
for relatively small values of c (see the curves for 
a = 7.3 ■ 10~^ and a = 6.6). By decreasing the 
value of a, the relaxation time increases consid- 
erably and the trapping of the particle occurs for 
lower values of c (see the curve for a = 2.7-10-2). 
For the lowest value, a = 10"^, we recover the 
regime of strong fluctuations due to the predomi- 
nance of the noise term with respect to the revert 
term in Eq. (8). The fluctuating behavior of all 
the curves before the divergence in Fig. [5] is also 
due to this effect. 

It is interesting to show, for our NLH model 
(Eqs. (7)-(9)), some of the well-established sta- 
tistical features of the financial time series, such 
as the probability function (PF) of the stock 
price returns, the PF of the volatility and the re- 
turn correlation, and to compare them with the 
same characteristics obtained from real market 
data. As initial position we choose Xq = —0.75. 
For this start point, located inside the well, we 
have very interesting behavior of the MET as 
a function of the noise intensity b (see panel 
(a) of Fig SIA). In Fig. [6] we show the PF of 
the returns. To characterize quantitatively this 
PF with regard to the width, the asymmetry 
and the fatness of the distribution, we calcu- 
late the mean value < Ax >, the variance 
a Ax, the skewness K3, and the kurtosis K4 for 



the NLH model and for real market data. We 
obtain for real data: < Ax >= 0.221 ■ 10"^, 
(TAx = 0.0221, K3 = -1.353, K4 = 79.334; 
for the NLH model: < Ax >= -1.909 • 10"^ 
(TAx = 0.0246, K3 = -4.3501, K4 = 441.65. The 
agreement between theoretical results and real 
data is quite good except at high values of the 
returns. These statistical quantities clearly show 
the asymmetry of the distribution and its lep- 
tokurtic nature observed in real market data. In 
fact, the empirical PF is characterized by a nar- 
row and large maximum, and fat tails in compar- 
ison with the Gaussian distribution [Mantegna 
l 




FIG. 6 Probability function of the stock price 
returns: (a) real data (circle), (b) NLH model 
(Eqs. (7)-(9)) (triangle). The values of the param- 
eters are: a = 2.00, b = 0.01, c = 0.75, a;o = -0.75, 

Xabs = -6.0, Vstart = 8.62 X 10"^. 

Fig. [7] we show the PF of the volatility for our 
model, and we can see a log-normal behavior 
as that observed approximately in real market 
data. The agreement is very good for values of 
the volatility greater than v ~ 0.03 and discrep- 
ancies are in the range of very low volatility val- 
ues. Finally in Fig. [8^ we show the correlation 
function of the returns for NLH model and for 
the real data. The agreement of the two corre- 
lation functions is very good for all the time. As 
we can see, the autocorrelations of the asset re- 
turns are insignificant, except for very small time 
scale for which microstructure effects come into 
play. This is in agreement with one of the styl- 
ized empirical facts emerging from the statistical 
analysis of price variations in various types of fi- 
nancial markets [Cont, 2001]. The good agree- 
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FIG. 7 Probability function of the volatility: (a) real 
data (black circle), (b) NLH model (Eqs. (7)-(9)) 
(solid line). The values of the parameters are the 
same as in Fig. [51 



ment between theoretical and experimental be- 
haviour is confirmed by the correlation function 
of the logarithmic absolute returns, which de- 
cays slowly to zero (see Fig. [8)3). Our last in- 
vestigation concerns the PF of the escape time 
of the returns, which is the main focus of our 
paper. By using our model (Eqs. (7)-(9)), we 
calculate the probability function for the escape 
time of the returns. We define two thresholds, 
Axi and Ax/, which represent the start point 
and the end point respectively for calculating 
MET. To fix the values of the two thresholds 
we consider the standard deviation (SD) a^x of 
the return series over a long time period cor- 
responding to that of the real data and we set 
Axi = — O.lcJAa;! Ax/ = —1.5 a Ax- The initial 
position is xq = —0.75 and the absorbing bar- 
rier is at Xabs = —6.0. We use a trial and er- 
ror procedure to select the values of the param- 
eters a, b, and c for which we obtain the best 
fitting between all the statistical features con- 
sidered, theoretical (NLH model) and empirical 
(real data). As real data we use the daily clo- 
sure prices for 1071 stocks traded at the NYSE 
and continuously present in the 12— year period 
1987 - 1998 (3030 trading days). In Fig. M we 
report the results for the PF of the escape times 
obtained both from real and theoretical data: 
we note a good qualitative agreement between 
the two PFs. Moreover we check the agreement 
between the two data sets by performing both 
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FIG. 8 (a) Correlation function of returns (panel a) 
and log absolute returns (panel b): real data (circle), 
NLH model (Eqs. ^ - ^) (solid line). Inset: detail 
of the behaviour at short times. The values of the 
parameters are the same of Fig. [51 Inset: detail of 
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FIG. 9 Probability function of the escape time of the 
returns from simulation (solid line), and from real 
data (black circle). The values of the thresholds are: 
Axi ~ — O.Ictak, Axf = — l.ScTAa;- The values of the 
other parameters are the same as in Fig. [S] 
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and Kolmogorov-Smirnov (K-S) goodness- 
of-fit tests. The results are: = 0.01620, 

= 0.00017 (x^ indicates the reduced x^) and 
D = 0.14, P = 0.261, where D and P are re- 
spectively the maximum difference between the 
cumulative distributions and the corresponding 
probability for the K-S test. The results ob- 
tained from both tests indicate that the two dis- 
tributions of Fig. [9] are not significantly different. 
Of course, a better quantitative fitting procedure 
could be done by considering also the potential 
parameters. This detailed analysis will be done 
in a forthcoming paper. 



IV. CONCLUSIONS 

We studied the statistical properties of the 
escape times in different models for stock mar- 
ket evolution. We compared the PFs of the es- 
cape times of the returns obtained by the basic 
geometric Brownian motion model and by two 
commonly used SV models (GARCH and Heston 
models) with the PF of real market data. Our 
results indeed show that to fit well the escape 
time distribution obtained from market data, it 
is necessary to take into account the stochas- 
ticity of the volatility. In the nonlinear Heston 
model, recently introduced by the authors, af- 
ter reviewing the role of the CIR parameters on 
the dynamics of the model, we analyze in detail 
the role of the initial conditions on the escape 
time from a metastable state. We found that 
the NFS effect, which could be considered as a 
measure of the stabilizing effect of the noise on 
the marked dynamics, is more pronounced for 
unstable initial positions near the maximum of 
the potential. For initial positions inside the po- 
tential well we recover an interesting nonmono- 
tonic behavior with a minimum and a maximum 
for the MET as a function of the parameter b. 
This behaviour is a typical signature of the NFS 
effect in the divergent dynamical regime [Fias- 
conaro et al, 2005]. To check the reliability of 
our NLH model we compare the return correla- 
tion function, the PFs of the returns, the volatil- 
ity and the escape times with the corresponding 
ones obtained from real market data. We find 
good agreement for some of these characteris- 



tics. 
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